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Abstract. We show that contrary to what is expected from ID stationary model atmospheres, 
3D hydrodynamical modeling predicts a considerable influence of convection on the spectral 
properties of late-type giants. This is due to the fact that convection overshoots into the for¬ 
mally stable outer atmospheric layers producing a notable granulation pattern in the 3D hydro- 
dynamical models, which has a direct influence on the observable spectra and colors. Within 
the framework of standard ID model atmospheres the average thermal stratification of the 3D 
hydro model can not be reproduced with any reasonable choice of the mixing length parameter 
and formulation of the turbulent pressure. The differences in individual photometric colors - in 
terms of 3D versus ID - reach up to ~ 0.2 mag, or ATeff ~ 70 K. We discuss the impact of full 
3D hydrodynamical models on the interpretation of observable properties of late-type giants, 
briefly mentioning problems and challenges which need to be solved for bringing these models 
to a routine use within the astronomical community in 5-10 years from now. 
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1. Introduction 

Convection plays an important role in governing the interior structure and evolution 
of late-type giants (i.e., stars on the red and asymptotic giant branches, RGB/AGB). 
Besides of aiding the energy transport from stellar interiors to the outer layers, convection 
is important in delivering heavy elements from the nuclear burning layers to the outer 
atmosphere. Since convective mixing changes the chemical composition both in the outer 
atmosphere and stellar interiors, it eventually alters also the stellar structure because of 
changes in atomic and molecular opacities. This affects observable properties of a star, 
lifetimes in different evolutionary stages, and so forth. Obviously, proper understanding 
of convection is of fundamental importance for building realistic evolutionary models, 
which form the basis of our understanding of individual stars and stellar populations. 

Convection in current theoretical models is treated in a rather simplistic way, typically 
within a framework of mixing length theory (MLT). Inevitably, this has a number of 
drawbacks. For instance, the efficiency of convection within the framework of MLT is 
scaled by a-priori unknown mixing-length parameter, Q!mlt (defined as the ratio of the 
mixing length to the pressure scale height), usually calibrated with Solar models. Almost 
certainly, ajuLT needs not to be the same in main-sequence stars, subgiants, giants and 
supergiants, as is assumed in current evolutionary models (see, e.g., Freytag & Salaris 
1999 for more details on this issue). Not surprisingly, there have been many attempts 
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Figure 1. Snapshot of the emergent white light intensity during the temporal evolution of a 
hydrodynamical red giant model. Note the spatial scales of the granulation pattern. The relative 
rms intensity contrast of the granulation pattern is 22.5 % at this particular instant in time. 

during the last few decades to improve the treatment of convection in stellar models (for 
example, implementing the concept of convective overshooting, which allows for convec¬ 
tion to penetrate beyond the classical boundaries of a convective layer). While such efforts 
are incremental steps towards a more realistic modeling of convection in stellar interiors, 
a fundamental breakthrough in this area is likely to go beyond the classical stationary 
ID modeling. Obviously, this may be possible with full 3D hydrodynamical models, as 
they account for time-dependent and three-dimensional character of convection from first 
principles, providing a degree of realism in the treatment of non-stationary phenomena 
(and convection in particular) that is beyond reach with classical approach. 


2. Convection in late-type giants with 3D hydrodynamical models 

To investigate the influence of convection on the observable properties of late-type 
giants we employed a radiation-hydrodynamics code named CD^BOLD (mainly developed 
by B. Freytag and M. Steffen; for a description of the code see Wedemeyer et al. 2004) 
to construct a fully time-dependent 3D hydrodynamical model of a prototypical late- 
type giant (Tefi ~ 3700K, \ogg = 1.0, [M/H] = 0.0; Fig. ^|. The hnal goal of this work 
was to make a comparison of the thermal structures and spectral properties of the 3D 
model with those predicted by classical ID model atmospheres. Note that ID atmosphere 
model used in this comparison was calculated using the same atmospheric parameters 
and employed the same physical input data (opacities, equation of state, description of 
radiative transfer) as the hydrodynamical model. 

Indeed, the thermal structure of the 3D model is quite different from that obtained with 
the classical ID approach (Fig.[21). Three ID models shown in Fig.|21were calculated for 
different mixing length parameters (omlt = 1-0,1.5,2.0) assuming a vanishing turbulent 
pressure, i.e. Pturb = 0. Two further ID models were constructed including a non-zero 
turbulent pressure, calculated as Pturb = fpv'^, where p and v are gas density and velocity 
respectively, and / is a dimensionless factor (usually / < 1). Interestingly, no reasonable 
combination of omlt and / allows to construct a classical ID model which reproduces 
the structure of 3D hydrodynamical model (Fig. [2). In the framework of MLT convection 
is confined to optically thick regions (logP > 3.5), and consequently the optically thin 







Convection in late-type giants 


3 



log,„P^.,[dyncm'"] 

Figure 2. Temperature stratification in the 3D hydrodynamical model of a late-type giant, as 
a fnnction of gas pressure, -Pgas- Thick solid line is a 3D hydrodynamical model (averaged on 
surfaces of equal geometrical depth), thin lines are ID plane-parallel models calculated using 
different mixing-length parameter (aniLT = TO, 1.5, 2.0, from left to right; in all cases Pturb = 0). 
Two ID models with non-zero turbulent pressure are given by dashed (omlt = 2.0 and / = 1.0) 
and dashed-dotted (oivilt = 2.0 and / = 2.0) lines. 

regions are hardly sensitive to convection. However, the geometric distance between the 
upper boundary of the convective region and optical depth unity is not large, and in 
the 3D hydrodynamical model convection overshoots into the optically thin layer. This 
produces a prominent granulation pattern on the stellar surface (cf. Fig. ^ which has a 
direct influence on the emerging spectrum, and thus - the photometric colors. 

The differences between colors predicted with 3D and ID models are significant, e.g., 
A{V - K) ^ 0.2, or AT^s ~ 70K (Fig. El see Kucinskas et al. 2005 for details on 
the calculation of spectra and colors with the 3D model). Note that this difference is 
comparable to a typical error margin in the individual Tefj-color relations (Kucinskas 
et al. 2006) which are used to transform theoretical isochrones to color-magnitude planes. 


3. Conclusions 

Clearly, our 3D model of a prototypical late-type giant predicts considerably different 
thermal structures from those inferred with the classical ID model atmospheres. Spec¬ 
tral properties of the 3D model are rather sensitive to convection too, which results in 
significant differences between photometric colors calculated with the 3D and ID model 
atmospheres. All this should be properly taken into account both with evolutionary mod¬ 
els and model atmospheres of late-type giants. At the same time, there is a number of 
issues that will have to be tackled to improve the 3D models, just to name a few: 

• Implement a possibility of direct spectral synthesis with 3D stellar atmosphere mod¬ 
els, preferably in non-LTE; 

• Improve spectral line databases, atomic and (especially!) molecular opacities; 

• Investigate the properties of stars in the regions of HR diagram that were poorly 
covered with 3D models up to now, at different metallicities (especially RGB/AGB stars). 

While 3D hydrodynamical models are still computationally expensive today, grids of 
synthetic spectra calculated in 3D approach may be available in 5-10 years from now (see 
Ludwig & Kucinskas 2005 for a discussion). No doubt, this will open new possibilities 
for improving models of stellar evolution (especially those on RGB/AGB), and will alow 
to study a variety of new phenomena that are beyond reach with classical ID models. 
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Figure 3. Influence of surface granulation on the broad-band photometric colors of red giant, 
as reflected by magnitude differences in various band-passes (indicated on the top of the panel) 
between the predictions of 3D hydrodynamical and classical ID model atmospheres. The different 
lines depict various approximations employed in the treatment of the 3D radiative transfer (see 
Kucinskas et at 2005 for details). 
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Discussion 

Feast: At least in case of large amplitude variables one can circumvent the models 
by empirical calibration of luminosity, age and metallicity in our Galaxy and the Local 
Group. 

Kucinskas: Galibration of selected stellar properties is indeed possible with observa¬ 
tions. However, all three quantities you mention are derived employing certain theoretical 
assumptions, thus they are not entirely model-independent. 

Feast: What is it in your models which determines the size of the convective cells? 

Kucinskas: To first approximation the size of the convective cells is set by (five to ten 
times) the pressure scale height at the stellar surface. In earlier spectral types with more 
intense radiation fields the cell size is somewhat biased towards larger values. Sphericity 
related effects in late-type supergiants are not fully clarified yet. However, there exist 
indications of a bias again towards larger cells. 

Ardeberg: With the new synthetic spectra and colors available, should we foresee any 
problems related to our current models, simulations and conclusions? 

Kucinskas: I would expect that predictions of the current theoretical models will not 
change drastically. There is no doubt, however, that certain details of the global picture 
will be different. How much - this needs to be clarified with further modeling. 




